Some Integrable Quantum Systems on the Lattice 
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Abstract 

The Weyl relations, the harmonic oscillator, the hydrogen atom, the Dirac equa- 
tion on the lattice are presented with the help of the difference equations and the 
orthogonal polynomials of discrete variable. This area of research is attracting more 
interest due to the lattice field theories and the hypothesis of a finite space. 



1 Weyl group on finite space 

We defined the position space of dimension N with orthonormal basis 



|0) 



/1\ 





1 





\N-1) 








I j) = k 



Similarly we construct an N-dimensional momentum space with ortonormal basis 
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where uo = 1, and (/ | A;) = 5/ 
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Two operators acting on these spaces are defined as 
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from which we construct Ua = A"-, Vi, = B^, a,b E Z 

On the position space, we have Ua \j) — \j — a) , H \j) — lo^^ (mod. N) and in the 
momentum space Ua\k) — uj^^ j/c) , Vh\k) — \k + h) , (mod. N) 

We can define the representation of vectors and operators as follows. 

From the expansion in the momentum space 

1 N-l N-l 

\F) = J^Y.\k) {k\F)=Y.h,\k) 

fe=0 fc=0 

we obtain the representation in position space 

F{j)^{3\F)^Y.^,h{3) (1) 

where/fc (j) = (j | k) = 

Similarly, from the expansion in momentum space. 
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we obtain the representation in momentum space 

G{k)^{k\G)=Y.ajf,{k) 
For the operators T4, Vf, we have in position space 



(2) 



(UaF) {j) = {j\Ua\F) = Fij + a) 
iV,F) {j)^{j\V,\F)^u;-'^F{j) 



and in momentum space 



(UaG) (k) = {k\Ua\G)=UJ-''''G{k) 

(HG) {k) = {k\V,\G) = G{k + b) 
From (1) and (2) we construct a Finite Fourier transform 
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The Weyl approach to Quantum Mechanics[l] is based in the properties of the opera- 
tors A, B when the N-dimensional space becomes infinite. 

Postulate I. There exist two parameter abehan group in an N-dimensional space whose 
elements A and B satisfy 

AB = uBA , u^' = 1 
A^ \j) = \j + s) B' \3) = u^' \3) 
Postulate II. In the continuous limit A?" — > oo we can identify 

uj = e'^ — > e*^'', ^ << 1, 77 << 1 

A' = (^e'^^y — >e'''^, ^s^a 

B' = (e^"^)* e'^^, rjt^T 

A'B* = u'^B^A' — > e'^^e'^^ = e'^^e'^^e'"^ 

The justification of Postulate II lyes in the fact that in the continuous case the action 
of the translation and multiplication operators is 

U^f(x)^e"'''fix) = fix + a) 

VJix)^e^^Qfix) = e'^''f{x) 

which are equivalent to the relations of Postulate I. From these equations the interpreta- 
tion of the operators P, Q is derived as the generators of the translations and multiplica- 
tions operators. 



2 The harmonic oscillator on the lattice 

In the discrete case we take the Kravchuk polynomials k^^^ {x, N) with a; = 0, 1, 2, ■ ■ • 

— 1, [2] and for the normalized functions the Wigner functions, that appear in the 
representation of the rotation group, dl^^, 

with N = 2j m = j — n m' = j — x p = sin^ — , Q = cos^ — 

After substitution in the fundamental formulas for the orthogonal polinomials we get 
[3] for the creation and annihilation operators 

^p(N -x-n) d]_^^ {13) + ^Jpqx {N - x + l)o?j_„_ 
V^p?n(A^-n + l)4_„+i, {(3) = 

= p{N-x-n) d]_^^ {(3) + ^pq{x + l) {N - x)d]_^^ {(3) 



The last equations can be written down in terms of the new parameters 

i sin (5^ {j + m){j-m+ l)dl_^^^, = sin^f (m + m') d^^^, + 

+isin - m') {j + m' + l)dl^^,^, {(5) 



\ sin N {j -m){j + m + lX+i,m' (/?) = sin^l + m') dLm' iP) + 



+|sin + m') {j -m' + l)<^,_i 



The creation and annihilation operators are connected with the raising and lowering 
operators for the spherical harmonics Yjm- In fact, from the connection between rf^^/ and 
Yjm, we get 

^^^"^ ^ - m) (j + m + = ;^'^+^^'"^ 



^^-^^■"^ " + (j + ^ + l)>S-,m-l = ^>^-^Jm 

In order to make more transparent the connection between the creation and annihi- 
lation operators with the raising and lowering operators of the spherical harmonics, we 
take the commutation and anticommutation relations of the former operators. 

wich in the limit j — > oo goes to 

a, al'] -^n (s) = V'n (s) 

Similarly 

(^^t + a^a) O = |(2n + 1) - ^1 d?^^, {(3) 
which in the limit j — > oo goes to 

(aa^ + a^aj il^n (s) = (2n + 1) ipn (s) 

If we multiply both sides by huj/2 we obtain the eigenvalue equation for the hamilto- 
nian. 

The interpretation of this model can be taken from the quantum harmonic oscillator. 

The energy levels are equally distant by the amount hur and are labelled by n = 
0, 1, 2, • • • oo. In the quantum harmonic oscillator of discrete variable we have also the 
discrete eigenvalues of the hamiltonian connected with the index m — j — n oi the Wigner 
function d^^, 

These values are equally separeted but finite (m = —j, ■ ■ ■ + j). Similarly the eigenval- 
ues of the position operator A + A~^ are also discrete and connected to the index m' = j —x 
of the Wigner functions but finite {m' — —j, ■ ■ ■ 



The integer numbers ,t = 0, 1, • ■ ■ 2j are related to the quantity x = as where s is the 
continuous variable and a = ^Mu/h. Since a; is a pure number and s has the dimension of 

a length, the spacing of the one-dimensional lattice is equal to 1/a = ^Jh/Muj. Therefore 
the Planck's constant h play role with respect to discrete space similar to the role with 
respect to discrete energy values. 



3 The Hidrogen atom in the lattice 

We start from the difference equation for the Meixner polynomials, the limit of which 
goes to the Laguerre polynomials in the continuous case. For the normalized Meixner 
polynomilas M^'^ (x) = d~^Jp{x)m^ [x] we get the following difference equation [3] 



^J^I{^i^^x)~{^+l)M,, (a; + 1) + ^px (x + 7 - 1)M„ [x-l)- 
- [// (x + n + 7) - n + x] M„ (x) = 

For the Hidrogen atom in the continuous case one takes as the solution of the reduced 
radical equation the wave function 

is) = yfp^)Ll^U is) , Pi (s) = sp (s) 
where p (s) is the weight function. In the discrete case we take the wave function 

Un(x) ^ d~^^/p^M2{x) , pi^ p(x + -f)p(x) 
The difference equation now reads: 



J^^^U„ (X + 1) + ,/^f/„ (X - 1) - t^fE±jl±lu„ (.) = - 1) 

Vx-l-7-l-l yx-1-7 X-I-7 X + J 

This equation is of Sturn-Liouville type, from which an orthogonality relation can be 
derived: 

^ Urn jx) Un jx) „ , 

2^ , = 0, if n / m 

We can construct also raising and lowering operators for the normalized Meixner 
functions. We get 

L+C4 (x) = ^/i (7 + n) {n - l)[/„+i (x) = (a; + n + 7) [/„ (x) - px {x + -f)Un {x - 1) 



I / X -\- \ 

L~Un (x) = Jpn{n + 'y - l)Un-i = p{x + 'y + n)Un {x)-p {x + 7) W — ■ —Un {x + 1) 

* V ^ + 7 + 

The action of these operators is to create or annihilate a new state the eigenvalue of 
which (with respect to the energy operator) is increased or decreased by unity. 



4 Dirac and Klein-Gordon equation on the lattice 

Prom the Dirac equation in momentum space {'j/iPfi — itlqc) ip {p) =0 we can construct 
the wave equation in position space with the help of the Fourier transform. We define the 
following difference operators 

A/ (j) = / (j + 1) - / (j) , A/ (j) = 1 {/ 0- + 1) + / (j)} , 

v/(j) = /(j)-/(j-i), v/(j) = ^{/(j) + /(j-i)}, 

and the partial difference operators with respect to a function of several discrete variables 

Az./ {j,,) = / {j^ + S,,^) - f {j^) , 

and similarly V^f (j^) and V^f (j/x)- 
Prom these operators we construct 

^ Ja, n A. , s; ^ iv, n V. , 

3 3 

^ n A, , ^- ^ n V, . 

From the Fourier transform we can derive the wave equation in lattice space. 

The kernel of the transform satisfies: 

1 2 ~ 

- exp (27ri (/c.j) e) = i- tan {irk^e) exp (27ri (k.j) e) 

We could apply the Fourier transform to the Dirac equation in momentum space and 
we would obtain the discrete wave equation. Instead we postulate a difference equation 
that in the limit goes to the continuous differential equation, namely, 

{^Y^t^ - mocri+^ ip {jf,) = 0. 

The kernel of the Fourier transform or "plane wave" is a particular solution of this 
equation if it satisfies 

^7'*- tanvr/c^e — rnocj r]~^ exp 27ri {k ■ j)s = 0. 

Applying the operator 7'*^ tan irk/^e + rriQC from the left to the last equation we obtain 

4 

— (tan nk^e) (tan vr/c^e) — m^c^ = 0, 

which is the intcgrability condition for the solution of the wave equation. 

Applying to the wave equation the operator i^j'^Sj^ + mQcrj~ from the left we obtain 
the discrete version of the Klein-Gordon equation in the lattice space 

a particular solution of which is again the "plane wave" provided the intcgrability condi- 
tions is satisfied. 



Acknowledgment (s) 



The autor expresses his gratitude to Prof. Doebner, Dobrev, Hennig and Luecke for the 
invitations to the Symposium. 

This work has been partially supported by D.G.I.C.Y.T. (Spain) ^ Pb 96 0538. 

References 

[1] H. Weyl, The Theory of Groups and Quantum Mechanics, Dover, New York 1950. 

[2] A.F. Nikiforov, S.K. Suslov, V.B. Uvarov, Classical Orthogonal Polynomials of a 
Discrete Variable, Springer, Berlin 1990. 

[3] M. Lorente, "Creation and Annihilation operators for orthogonal Polynomials of 
continuous and Discrete Variables", Electronic Transactions on Numerical Analysis, 
Vol. 9 (1999), pp. 102-111. 



[4] M. Lorente, "Representations of the Discrete Inhomogeneous Lorentz Group and 
Dirac Wave equation on the lattice", J.Phys.A: Math. Gen. 32 (1999) 2.481-2.497. 



